home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / sgecon.z / sgecon
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSGGGGEEEECCCCOOOONNNN((((3333FFFF))))                                                          SSSSGGGGEEEECCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SGECON - estimate the reciprocal of the condition number of a general
  10.      real matrix A, in either the 1-norm or the infinity-norm, using the LU
  11.      factorization computed by SGETRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE SGECON( NORM, N, A, LDA, ANORM, RCOND, WORK, IWORK, INFO )
  15.  
  16.          CHARACTER      NORM
  17.  
  18.          INTEGER        INFO, LDA, N
  19.  
  20.          REAL           ANORM, RCOND
  21.  
  22.          INTEGER        IWORK( * )
  23.  
  24.          REAL           A( LDA, * ), WORK( * )
  25.  
  26. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  27.      SGECON estimates the reciprocal of the condition number of a general real
  28.      matrix A, in either the 1-norm or the infinity-norm, using the LU
  29.      factorization computed by SGETRF.
  30.  
  31.      An estimate is obtained for norm(inv(A)), and the reciprocal of the
  32.      condition number is computed as
  33.         RCOND = 1 / ( norm(A) * norm(inv(A)) ).
  34.  
  35.  
  36. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  37.      NORM    (input) CHARACTER*1
  38.              Specifies whether the 1-norm condition number or the infinity-
  39.              norm condition number is required:
  40.              = '1' or 'O':  1-norm;
  41.              = 'I':         Infinity-norm.
  42.  
  43.      N       (input) INTEGER
  44.              The order of the matrix A.  N >= 0.
  45.  
  46.      A       (input) REAL array, dimension (LDA,N)
  47.              The factors L and U from the factorization A = P*L*U as computed
  48.              by SGETRF.
  49.  
  50.      LDA     (input) INTEGER
  51.              The leading dimension of the array A.  LDA >= max(1,N).
  52.  
  53.      ANORM   (input) REAL
  54.              If NORM = '1' or 'O', the 1-norm of the original matrix A.  If
  55.              NORM = 'I', the infinity-norm of the original matrix A.
  56.  
  57.      RCOND   (output) REAL
  58.              The reciprocal of the condition number of the matrix A, computed
  59.              as RCOND = 1/(norm(A) * norm(inv(A))).
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSGGGGEEEECCCCOOOONNNN((((3333FFFF))))                                                          SSSSGGGGEEEECCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      WORK    (workspace) REAL array, dimension (4*N)
  75.  
  76.      IWORK   (workspace) INTEGER array, dimension (N)
  77.  
  78.      INFO    (output) INTEGER
  79.              = 0:  successful exit
  80.              < 0:  if INFO = -i, the i-th argument had an illegal value
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.